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ABSTRACT 

We present high spatial resolution mid- and far-infrared images of the Vega debris disk obtained 
with the Multiband Imaging Photometer for Spitzer (MIPS). The disk is well resolved and its angular 
size is much larger than found previously. The radius of the disk is at least 43" (330 AU), 70" (543 
AU), and 105" (815 AU) in extent at 24, 70 and 160 /im, respectively. The disk images are circular, 
smooth and without dumpiness at all three wavelengths. The radial surface brightness profiles follow 
radial power laws of r -3 or r~ 4 , and imply an inner boundary at a radius of 11"±2" (86 AU). 
Assuming an amalgam of amorphous silicate and carbonaceous grains, the disk can be modeled as 
an axially symmetric and geometrically thin disk, viewed face-on, with the surface particle number 
density following an inverse radial power law. The disk radiometric properties are consistent with a 
range of models using grains of sizes ~1 to ~50 /im. The exact minimum and maximum grain size 
limits depend on the adopted grain composition. However, all these models require a r _1 surface 
number density profile and a total mass of 3±1.5xl0 -3 Mq in grains. We find that a ring, containing 
grains larger than 180 /im and at radii of 86-200 AU from the star, can reproduce the observed 850 
/im flux, while its emission does not violate the observed MIPS profiles. This ring could be associated 
with a population of larger asteroidal bodies analogous to our own Kuiper Belt. Cascades of collisions 
starting with encounters among these large bodies in the ring produce the small debris that is blown 
outward by radiation pressure to much larger distances where we detect its thermal emission. The 
relatively short lifetime (< 1000 years) of these small grains and the observed total mass, ~3xl0 -3 
M®, set a lower limit on the dust production rate, ~10 15 g/s. This rate would require a very massive 
asteroidal reservoir for the dust to be produced in a steady state throughout Vega's life. Instead, we 
suggest that the disk we imaged is ephemeral and that we are witnessing the aftermath of a large and 
relatively recent collisional event, and subsequent collisional cascade. 

Subject headings: circumstellar matter - infrared: stars - planetary systems - stars: individual (Vega) 



1. INTRODUCTION 

Vega (a Lyrae = HD 172167 = HR 7001, AO V, d=7.76 
pc) has been used both as a fundamental photometric 
standard and as a template for modeling of stellar at- 
mospheres. One of the highlights of the IRAS mission 
was the discovery of a large infrared excess around Vega, 
which was observed originally as a photometric stan- 
dard. The far infrared emission was attributed to ther- 
mal dust emission fro m a disk of debris orbiting the star 
ijAumann et al.l 119841) . Subsequent IRAS observations 
of other main sequence stars in the solar neighborhood 
found that nearly 15% of them exhibit infrared exces ses 
ijBackman fc ParescdlT993t Mannings fc BarlowllT998|) . 

The dust found around these main sequence stars 
cannot be primordial material left over from the star 
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forming stage because the time scale to remove primor- 
dial material (<10 Myr) is short compared to the life- 
times of these stars (for example, Vega is estimated 
as 35 Myr old, iBarrado v Navascueslll998l iSong et all 
2000). Therefore, the dust must be re-supplied. Second- 
generation dust in such systems is thought to arise pri- 
marily from collisions between planetesimals and from 
cometary activity. Ground- and space-based observa- 
tions of the nearest subsample of Vega-like stars from 
optical to submillimeter wavelengths have revealed disk- 
like structures — debris disks — with density gaps and 
enhancements i jpc^Ila^c} el; alJ|1998l: |Gr£avgf^^djEj93 
[ Koerner. Sargent fc Ostroffl 120011: iWilner et alJ l2002t 
iWeinberger et alJ l2002t iClampin et alJ 12003^ Debris 
disks are the most visible signposts of other planetary 
systems, representing indirect evidence of planetary sys- 
tem formation. We can learn about the diversity of plan- 
etary systems from the study of debris disk structures. 
Vega is one of the closer and brighter debris disks, which 
provides us an opportunity to study it in detail spatially, 
as a foundation for unders t andin g other debris disks. 

Harvcv. Wilk ing. & Jo v (1984) observed Vega at 47 
and 95 /im with the Kuiper Airborne Observatory and in- 
ferred that the s ource of infrared emission may be as large 
as 23" in radius. iVan der Bliek. Prusti. fc Waters! (1994) 
re-analyzed the IRAS data and found that the IRAS 60 
/xm emission from Vega is 18"±3" in radius. They also 
suggested that the emission at 60 /im is caused by small 
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dust grains with size between 0.1 and 10 /im. Infrared 
Space Observatory (ISO) observations show a smooth, 
resolved, face-on disk wit h a radius of 22"±2" at 60 /im, 
and 36" ±3" at 90 um. (iHeinrichsen. Walker, fc Klaasl 
1 1998ft . iMauron fc Dole! (|1998l) attempted to detect 
Vega's disk through optical scattered light by using lin- 
ear photopolarimetry. Their upper limit implies that the 
debris around Ve ga contains a very smal l numb er of 0.01- 
0.3 /im grains. IZuckerman fc Becklinl l|1993j) detected 
emission by dust from Vega at 800 /jm and suggested 
the presence of 100-300 /im grains to account for the 
submillimeter observations. 

The 850 /im map of Vega obtained by iHolland et alJ 
( 1998) with the Submillimeter Common-User Bolome- 
ter Array (SCUBA) on the James Clerk Maxwell Tele- 
scope shows an extended, roughly circular structure with 
dimensions of 24"x21"(±2"). The 850 /mi resolved 
emission is not uniformly distributed; instead, there is 
an elongated bright central region oriented northeast- 
southwest with a bright peak offset ~9" (70 AU) to the 
northeast of the star's position. Observat i ons at 1.3 mm 
obtained bvlKoerner. Sargent, fc Ostroffl l)2001|) and by 
IWilner et alJ l|2002j) resolved dust emission peaks offset 
from the star by 8" to 14" that appear to be associated 
with a ring of emission at a radius of 60 to 95 AU. 

Here we present the first high spatial resolution 24, 70 
and 160 (im ima ges of Vega obtaine d with the Spitzer 
Space Telescope ijWerner et al J 12004]) . The sensitivity 
and resolution of Spitzer allow us to detect fainter and 
more extended emission in the Vega debris disk system 
than was known previously. Details about observations, 
data reduction, and stellar photosphere subtraction are 
presented in The resolved images are discussed in 
The true nature of the disk can be revealed by studying 
the radial profiles of the disk, which are presented in 21 
with interpretations of disk temperature structures and 
masses in SJS1 The origin of the debris we detected with 
Spitzer is discussed in ^1 followed by our conclusions in 

m 

2. OBSERVATIONS, DATA REDUCTION, AND 
PHOTOSPHERIC SUBTRACTIONS 

Vega was imaged on 11 April 2004 with the Multi- 
band I maging Photometer for Spitzer (MlPS: lRieke et aD 
2004a). Multiple dithered 3 sec exposures were obtained 
using the 24, 70 (both coarse and fine scales), and 160 
/im channels. Details about the observations are listed 
in Table ^ along with the combined sensitivity at each 
band. 

The data were processed using the MIPS instru- 
ment te am Data Analysis Tool (DAT; Gordon et al. 
I2004albf) for basic reduction (dark subtraction, flat field- 
ing/illumination correction). In addition, the 24 /im im- 
ages were processed to remove a vertical "jail-bar" pat- 
tern caused by the bright saturated central star. Due 
to the known transient behaviors associated with the 
70 /im detectors, a detector dependent structure in the 
coarse scale mode data was removed by subtracting col- 
umn averages from each exposure with the source re- 
gion masked. Similar artifacts were taken out in the 70 
/im fine mode data by subtracting the off-source chopped 
background observations. The processed exposures were 
then mosaicked together. 

The Vega disk at 24 /tm is tenuous compared to the 



very bright photosphere (7.4 Jy). The existence of a 
very low surface brightness disk is best illustrated by ex- 
amining the radial profile (Figure ^| . Relative to the 
observed Point Spread Function (PSF) scaled to Vega's 
photospheric level, the observed radial profile lies consis- 
tently above the scaled point source profile. The contrast 
between the bright and dark Airy rings in the image of 
Vega is much lower, implying the brightness of the dark 
Airy rings has been filled in by a low surface brightness, 
extended component. 

Vega's photosphere contributes a significant fraction of 
the flux near the center of our images at both 24 and 70 
/im, while the emission from Vega's debris disk dominates 
most of the flux at 160 /tm. To study the emission from 
Vega's debris disk the photospheric contributions have to 
be removed in all three channels. To remove Vega's pho- 
tosphere, we subtracted reference observed PSFs, scaled 
to Vega's photospheric flux, from our images. Chang- 
ing the PSF scaling factor by more than 10% resulted in 
artifacts in the subtracted images at the position of the 
core of the PSF and/or at the positions of the Airy min- 
ima and maxima. We conclude that our photospheric 
removal is accurate to 10% within any given resolution 
element in the resulting photosphere-subtracted images. 

At 24 /im we processed the data for a PSF reference 
star (HD 217382) in exactly the same way as we did 
for Vega. Because the central pixel in the Vega im- 
ages was saturated, registering the reference PSF to the 
Vega image was difficult. We took advantage of image 
latents present in both data sets to do the registration 
and photosphere subtraction on a frame-by-frame basis. 
We then mosaicked the resulting images to produce the 
final 24 /im image of Vega's disk, shown in Figure [21 
The mosaic does not show any saturation artifact at the 
center because of the sub-pixel dithers used in the ob- 
servation: a pixel that was saturated in one image is 
partially overlapped by unsaturated neighboring pixels 
in an image from a different dither position. There were 
enough dither positions to provide partial coverage even 
at Vega's position in the final mosaic. Because of this in- 
filling of the saturated area, the effective exposure time 
in the central 6" is about half of that elsewhere. Photo- 
spheric removal at 70 /im was more straight forward. We 
registered the scaled reference PSF mosaic (HD 48915) 
to the Vega mosaic by ccntroiding, and subtract ed. The 
results are shown in Figure |3a] (coarse scale) and !3bl ffine 
scale) . 

The MIPS 160 /tm array suffers from a spectral leak 
caused by an internal reflection in the optical train al- 
lowing leakage from very blue and bright Rayleigh-Jeans 
sources to contaminate the signals at 160 /tm. However, 
the spectral leak image is offset to one side of the true 160 
/tm image, and the brightness of the leak is proportional 
to the photospheric flux. Comparison of 160 /tm images 
of stellar (blue) sources with images of asteroids shows 
that the near-infrared leak contributes very little to the 
160 /im images on the opposite side of the source lo- 
cation. The predicted flux for Vega's photosphere is 162 
mjy, which is much fainter than the expected disk bright- 
ness at 160 /tm. We took advantage of this situation by 
using only the half of the 160 /tm image where the leak 
contribution is negligible. We also subtracted a scaled 
(red) 160 /tm reference PSF (asteroid Harmonia) from 
the Vega mosaic, using the pointing information (accu- 
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Fig. 1. — Radial profile of Vega (photosphere + disk) at 24 (im af- 
ter background subtraction. For comparison, the radial profile of 
an observed point source is scaled to Vega's photosphere flux and 
shown as filled circles. The observed radial profile lies consistently 
above the point source profile. The contrast between the bright 
and dark Airy rings in the PSF signatures is not as prominent as 
the one in the point source, suggesting the existence of a tenuous 
disk. 



rate to < 1") to register. The result is shown in Figure 

El 

3. DISK MORPHOLOGIES AT 24, 70 AND 160 /jm 

We define our observed sensitivity based on the 1-a 
background noise per pixel using blank sky area in the 
image. The 1-a background noise in the PSF-subtracted 
image is 11 /iJy/arcsec 2 at 24 /*m. The disk at 24 /jm is 
symmetric and centered at the star position; no obvious 
asymmetry is seen in the image. At the 1-a level, the 24 
/im disk extends to ~43" (330 AU) in radius. The total 
flux density (within the 1-a contour) is ~ 1.5 Jy ±10% 8 . 

This flux density value is in agreement with the IRAS 
25 fim measurement. The quoted IRAS 25 /im flux den- 
sity for Vega is ~10.5 Jy (combining IRAS PSC and 
FSC). Based on Kurucz models, Vega's photosphere is 
—6.63 Jy at 25 /im. The relation between the IRAS 
quoted flux and actual flux is F ac t ua i = F qU oted/K, where 
K is the color correction factor. The color correction fac- 
tor is 1.41 for T = 9750 K blackbody (Vega photosphere) 
and 0.83 for T = 95 K (debris disk as determined by pre- 
vious studies), suggesting that the actual flux from the 
disk is ~ 1.4±10% Jy at 25 /mi. 

Vega's disk looks very similar in both 70 fim images. 
The fine scale data were taken in 12 dithered positions 
with sub-pixel offsets to provide a good spatial sampling 
for the disk. There is no obvious clumpy structure in the 
fine scale image. Similar to the disk image at 24 /im, the 
disk at 70 /mi appears symmetric and smooth. At a 1-a 
detection limit, the outermost boundary of the disk at 
70 /mi (coarse scale mode) is —70" (543 AU) in radius. 
The total flux (within the 1-a contour) is -7 Jy (±20%). 

Assuming, based on the 24 and 70 /mi morphology, 
that the disk is also azimuthally symmetric at 160 /im, 
we are able to determine that the outermost extent of 
the disk at 160 /im is -105" (815 AU) in radius (1-a). 
Assuming the bad side of the disk image (affected by 
the leak) has similar flux as the good side, the total flux 

8 10% including errors in absolute flux calibration, and in color 
correction (less than 5% for a blackbody temperature of 95 K). 




Fig. 2. — Vega disk at 24 fim displayed with logarithmic scaling. 
North is up and east is toward the left. Due to the saturation of the 
central star at 24 fim, all the negative values after PSF subtraction 
have been excluded in the final mosaic, resulting in a smooth image 
at the core region (r < 6"). The coverage (effective exposure time) 
near the center is ~half of that outside the saturation region. The 
instrument beam size (FWHM) of 6" at 24 fim is shown as a white 
circle in the lower-right corner. 



from the disk at 160 /im (within the 1-a contour) is —4 
Jy ±20%. 

Within the resolution errors, the disk image is circu- 
lar at all wavelengths, suggesting a face-on disk. This 
behavior is consistent with Vega being a pole-on star 
(Gulliver. Hill. & Adelman 1994). Based on previous 
studies (IRAS, ISO and SCUBA), the spectral energy 
distribution (SED) of the disk can be characterized as a 
blackbody of T = 95 K with a A -1 emissivity. The fluxes 
measured by MIPS match the previously observed SED 
within 15%. 

A striking immediate result from the MIPS images is 
the physical size of the Vega disk. Submillimctcr and mil- 
limeter observations show a low surface brightness ring 
with radius —11", much smaller than the emission seen in 
the MIPS bands. I n contrast, MIPS image s of the Foma- 
lhaut debris disk l)Stapelfeldt et all 120041) indicate that 
the submillimeter and far-infrared morphologies agree 
well. A second surprise is the presence of material warm 
enough to be detecte d in the MIPS bands at a large dis- 
tance from the star. lArtvmowicz &: Clampinl l)1997|) es- 
timated the size of the disk if the grains are blown out 
by radiation pressure, and suggested that the disk could 
be up to a few thousand AU, consistent with our images. 

4. THE SURFACE BRIGHTNESS RADIAL DISTRIBUTION 

4.1. At 24 /im 

Because the disk is almost perfectly face-on, its extent 
can be best shown in an azimuthally averaged radial in- 
tensity profile. The advantage of using radial profiles 
is that the noise is reduced by y/~N, where N is the to- 
tal number of pixels used in averaging the intensity at a 
given radius. The radial surface brightness profile of the 
disk at 24 /im is shown in Figure [5] A bin size of 2" 5 
(1 pixel) was used to calculate the average brightness at 
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Fig. 3a. — Vega disk at 70 (im in coarse scale mode displayed with 
logarithmic scaling. North is up and east is toward the left. The 
Vega photosphere has been subtracted off by scaling an observed 
PSF star. The instrument beam size (FWHM) of 18" at 70 ^m is 
shown as a white circle in the lower-right corner. 



a given radius r when r < 40", but a bin size of 10" (4 
pixels) was used when r > 40" to increase the signal-to- 
noise, up to 250" (~ 4') from the star. A power law plus 
a constant background (S(r) = S% 4 + r~ a ) was used to 
fit the data points between 22" and 150" to determine 
the true sky background (S% 4 ) and the flux radial depen- 
dence. A power index of — 4.1±0.1 is indicated, based on 
a least-squares fit. 

The dashed lines on Figure El mark the data range for 
the power law fitting. The power index of —4.1 fits fairly 
well for the disk surface brightness when r > 200 AU, 
but deviates significantly when r < 150 AU. To under- 
stand the inner structure of the disk, we first compute 
a disk model with a r~ a surface brightness distribution 
and then convolve the model disk with the 24 /im PSF. 
With a rapid power law fall-off (a =4), the radial profile 
of the convolved model disk (shown as a green solid line 
in FigureEJ is similar to the point source PSF, which was 
not seen in the observed data. This comparison suggests 
the surface brightness distribution must flatten or go to 
zero (a real physical empty hole) in the inner part of the 
disk. 

We first investigated the possibility of a real empty 
hole in the surface brightness distribution. By adjusting 
power-law indices and inner hole sizes, we are able to 
reproduce the observed radial profile using two different 
models. Model A, which is r~ 4 + S^ 4 with an inner disk 
hole (ri n ) at 5.5" with an error of ±7% at 90% % 2 con- 
fidence level (shown as a red solid line in Figure EJ), fits 
very well when r > 200 AU as expected. Model B, which 
is r~ 3 +S^ 4 with r in = 2.9" with an error of ±17% (shown 
as a blue solid line in Figure EJ), fits very well when r < 
200 AU, but has slightly higher surface brightness at r > 
200 AU compared to the r -4 power law. The lower panel 
of Figure El shows the ratio between the model and ob- 
served radial profiles. 

We have shown that a real empty hole in the surface 
brightness profile can fit the data well. Next we investi- 




FlG. 3b. — Same as Figure l3*al but in fine scale mode. 




Fig. 4. — Vega disk (half) at 160 /im displayed with logarithmic 
scaling. The right side of the disk is contaminated by the spectral 
leak, and therefore not shown in the image. The white asterisk 
at the center marks the stellar position. The instrument beam 
size (FWHM) of 40" at 160 (im is shown as a white circle in the 
lower-right corner. 



gated the possibility of a flat surface brightness distribu- 
tion in the inner part of the disk. Assuming the surface 
brightness is constant, S(ri n ) 7 when r < ri„, rather than 
zero in the empty hole case, we found two models that 
can fit the observed profile. Model A', which is r -4 + S^ 4 
with an ri n = 8.0" with an error of ±5% (shown as a red 
dashed line in Figure EJ , fits very well when r > 200 
AU. Model B', which is r" 3 + S^ 4 with r in = 4.4" with 
an error of ±14% (shown as a blue dashed line in Fig- 
ure El, fits very well when r < 200 AU, but has slightly 
higher surface brightness at r > 200 AU compared to the 
r~ 4 power law. The deficit in the flat surface brightness 
distribution is very large (>99%) when compared to the 
extrapolation of the steep power law to the disk center. 
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Fig. 5. — Radial profile of the Vega disk at 24 /im. A power law Fig. 7. — Radial profile of the Vega disk at 70 fim. A power law 

plus a constant background are used to fit the data points between plus a constant background are used to fit the data points between 

the two dotted lines. The best-fit power law index is —4.1 ± 0.1. the two dotted lines. The distribution is found to be consistent 

with a r~ 38 power law. 
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Fig. 6. — Radial profile of the Vega disk at 24 fim (open circles) 
compared with modeled surface brightness distributions after con- 
volution with the instrumental beam. Note that the background 
value determined in Figure|5]has been added back in the fits. The 
model fit with a steep power law (-4 or -3) and no inner boundary 
of the disk is shown as a green solid line, which resembles a PSF. 
The other models are: Model A — an empty hole with r~ 4 as a 
solid red line; Model B — an empty hole with r~ 3 as a solid blue 
line; Model A' - a flat distribution with r~ 4 as a dashed red line; 
Model B' — a flat distribution with r -3 as a dashed blue line. The 
lower panel of the plot shows the ratios between the modeled and 
the observed radial profiles. 



The difference between an empty hole and a flat distri- 
bution in terms of surface brightness is insignificant. For 
r > 200 AU, Model A and A' provide a good fit to the 
data whereas Model B and B' match the data well for 
r < 200 AU. The most important conclusion from these 
fittings is that the disk surface brightness distribution 
follows simple power-law dependences, implying that the 
disk density structure is simple and smooth. 

4.2. At 70 /jm 

The 70 /im radial profile of the Vega disk was computed 
using an average of the data from both observing modes 
for r <40", but only the coarse scale mode data when 
r > 40". The combined radial profile is shown in Figure 




Fig. 8. — Radial profile of the Vega disk at 70 fim compared 
with model surface brightness distributions after convolution with 
the instrumental beam. Similar to Figure|S] the background value 
has been added back in the model fits. Symbols and color scheme 
are the same as Figure |E| The lower panel of the plot shows the 
ratios between the modeled and the observed radial profiles. 



□ Fitting with a power law plus a constant background 
function, S(r) = S™ + r~", the surface brightness of 
the disk shows a r~ 3 - 8 power law for r > 200 AU with 
a very small error in the index. Similar approaches to 
those pursued at 24 /urn (hole or flat distribution) were 
used to fit the inner part of the radial profile at 70 /im. 
The results are shown in Figure |S1 For an empty hole in 
the surface brightness distribution, Model A gives a r,„ 
= 11.5" with an error of 5% and Model B give a r in = 
9.4" with an error of 7%. For a flat distribution, Model 
A' gives a ri n = 18" with an error of 6% and Model B' 
gives a r,„ = 16" with an error of 6%, and the deficit is 
larger than 99% when compared to the extrapolation of 
the power law to the disk center. For the inner part of 
the disk (r < 350 AU), Model A/ A' and Model B/B' give 
similar results, but Model A and A' give a much better 
fit for r >350 AU. 
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Fig. 9. — Radial dependent [24]— [70] distribution for the Vega 
disk. The observed color is plotted in filled squares. The color 
based on the canonical dust temperature is plotted as open dia- 
monds, compared to the ones derived from single sized silicates 
(a = 5.1 fim as the solid red line and a = 2 fim as the dashed red 
line). Th e [24]— [70] color derived from our two-component model 
(see 85.11 is also plotted as the blue dashed line. 



4.3. [24]- [70] color distribution 

To compare the disk structure between 24 and 70 fim, 
we first convolved the 24 /jm disk image with a Gaus- 
sion kernel to match the 70 /jm resolution. The radi- 
ally dependent [24]— [70] color is estimated at a given ra- 
dius from the star by computing the flux ratio between 
70 and 24 /im. The resultant radial color distribution 
is shown in Figure The [24]— [70] color initially in- 
creases (color temperature decreases) with distance from 
the star, but levels off between 300 and 600 AU (with 
a color temperature of ~67 K). For comparison, we also 
overplot the expected color from the canonical temper- 
ature distribution, T g = T*(R*/2rf- b with T* = 9750^ 
and i?» = 2.5i?0 for Vega, assuming spherical blackbody 
radiators (i.e., absorption efficiency Q a bs — 1)- 

We compute the radial color values for single-size sili- 
cates assuming the grains are heated only by stellar ra- 
diation and emit in an optically thin case (for details, 
see section EJ. As shown in Figure the [24]— [70] color 
in the inner part (r < 200 AU) of the disk is consistent 
with the emission from single-size silicates with radius 
a = 5.1 /im. Furthermore, a disk with only 2 /im grains 
does not reproduce the observed color distribution at all. 
The temperature structure of the disk cannot be eas- 
ily explained by thermal emission from a disk composed 
of single-size grains. To produce the observed [24]— [70] 
color, additional components must be invoked. 

4.4. At 160 /mi 

Since the disk is bright and extended at 160 /im, we 
were able to extract a disk radial profile based on the 
good half of the disk image. To verify our 160 /im radial 
profile is legitimate we also extracted radial profiles of a 
true (no leak) 160 fxxsx source (asteroid Harmonia) and a 
calibration star (HD 3712, KOIIIa) only using the good 
half of the image. We found that the radial profiles ex- 
tracted only from the good half the images look similar, 
suggesting the good half of the image is not affected by 
the leak. 

The Vega disk radial profile at 160 /im is shown in Fig- 
ure ^| after background subtraction. Similar approaches 
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Fig. 10. — Radial profile of the Vega disk at 160 /im compared 
with model surface brightness distributions after convolution with 
the instrumental beam. The background value is not included in 
the fits. Symbols and color scheme are the same as FigurelolandlSl 
The lower panel of the plot shows the ratios between the modeled 
and the observed radial profiles. 



(hole or flat distribution) were used to fit the radial pro- 
file. The results are: Model A (hole, r -4 ) gives a r,-„ 
= 15"; Model A' (flat, r" 4 ) gives a r in = 23"; Model B 
(hole, r~ 3 ) gives a r in = 14"; Model B' (flat, r~ 3 ) give 
a ri n = 21". The errors in these model fits are on the 
order of 30%. Again, the deficit in the flat distribution 
is greater than 99% when compared to the extrapolation 
of the steep power law to the disk center. 

4.5. Summary 

The surface brightness distribution at the 24 /im band 
agrees with a radial-dependent power law with index of 
—3 (inner part) or —4 (outer part). At 70 /im, the r 
power law fits better than the r in the outer part of 
the disk. In addition, this radial-dependent power law 
cannot apply all the way to the star; a hole or a flat 
distribution in the disk is required. The inner edge of 
the power-law surface brightness distribution depends on 
whether the central region contains no flux or a constant 
non-zero flux. In the latter case, the inner disk radius 
is ^50% larger. However, the current data are not ade- 
quate to distinguish between a empty hole and a flat non- 
zero central flux deficit. Therefore, the inner boundary 
of the disk is determined as a weighted average of these 
radii obtained from different fits in all three bands. We 
conclude that the inner boundary of the Vega debris disk 
has a radius of 11" ±2" (or 86±16 AU). The "hole" is 
not required to be physically empty; rather it is a region 
lacking material that emits strongly at 24, 70, and 160 
/im. Interestingly, the radius of the 11"±2" hole is about 
the same size as the 850 /im (and mm-wave) ring. This 
implies that the material we detect in the MIPS bands 
might originate from the submillimeter ring (large parent 
bodies) where collisional grinding generates small debris 
that is blown out by radiation pressure to the distances 
where we detect their emission. 
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5. DISK STRUCTURE AND MASS 

For particles moving in a steady-state flow, the mass 
gain/loss rate ^| through any annulus is expected to 
be constant. If the particles are moving at their termi- 
nal velocities (v ~ const.), the surface number density 
(£) is expected to be proportional to r _1 due to mass 
conservation (AM = 27rrwArE ~ const.). Therefore, for 
a radiation pressure driven outflow, the surface particle 
number density should follow an inverse radial power law 
if the particles have reached their terminal velocities. In 
H5.ll we will show the Vega debris disk is well modeled 
with an surface number density power law. 

5.1. Preferred Model 

We assume that the disk is axially symmetric and ge- 
ometrically thin, viewed face-on, and that the surface 
number density is governed by a simple radial power law; 
i.e., S(r) = Y,o(r/ro)~ p where r is the inner radius of 
the disk. Based on results in 21 we se t ?"o = 86 AU and 
simplify the disk inside ro as a physically empty space. 
The emission from grains at a given radius r and wave- 
length A can be written as 

2 

TTfl 

dF x = dn(r)-^-Q(a, \)B\(T r ), (1) 

where d is the distance to Vega and dn(r) — 27rS(r)rdr 
is the total number of grains of a given radius between r 
and r + dr. 

To probe the nature of the Vega disk further, we 
have fitted it with models including realistic grain op- 
tical properties. A dust grain temperature for a given 
set of grain parameters (composition, grain radius a, i.e. 
Qabs(a, A)) at a given wavelength (A) and location from 
the star (r), can be computed assuming the dust is only 
heated by the star. We use the photospheric flux from 
the Kurucz model for Vega as the input radiation field. 
The absorption coefficient is calculated with Mie the- 
ory, using the op tical constants for ast ronomical silicates 
(p g = 3.5 g/cm 3 . iLaor fc Draindll993|) . The grain tem- 
peratures at different radii from the star, T(r), shown 
in Figure 1111 are computed based on balancing the en- 
ergy between the absorption and emission by the dust 
grains (scattering can be ignored for the wavelengths of 
interest). 

The total flux from an annulus (r' to r") by integrating 
Equ. (1) is 

2 rr" 

F x (r') =2^o^ T Q(a,X)r p Q j f r^ +1 B x {T r )dr . (2) 

If the annulus is small enough (i.e., r' ~ r" ), the flux 
can be written as 

F x (r') = 2irX ^Q(a,\y (r')- p+1 B x (T r ,)(r" - r') . 

(3) 

The surface brightness distribution, S\(r'), is the total 
flux divided by the annulus area, i.e., 

S x (r') = 2S ^Q(a, A)rg ( r ')- p+1 Bx(T r ,)^^ . 

r (4) 

Equ. (4) can be used to qualitatively evaluate the 
expected power index in the surface number density. 
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Fig. 11. — Thermal equilibrium temperature distribution for 
different grains in the Vega environment. Astronomical silicates 
are plotted in blue, and amorphous carbon grains in red. Solid 
lines are for small grains (a = 2 fim) and dashed-dotted lines are 
for large grains (a = 18 /im). 



For a given wavelength and set of grain properties, the 
observed surface brightness, S(r), is proportional to 
B\(T r ) ■ r~ p . In general, the temperature, T r , is pro- 
portional to r -0 - 33 for small grains, or r -0 5 for large 
grains. Using the temperature distribution for silicate 
grains with a = 2 /im, the Planck function at 24 firn ba- 
sically follows r~ 2 in the inner part of the disk, but r~ 3 
in the outer part of the disk. In other words, if the ob- 
served 24 /im flux is dominated by the emission from 
small grains, a p = 1 power index for the surface density 
distribution is preferred based on the fits in 21 ( r_3 f° r 
inner disk but r~ 4 for outer disk). 

For proper comparison to the observed surface bright- 
nesses (after subtraction of the true sky background 
determined from the fits in 2J> the modeled surface 
brightness profile using Equ. (4) is then convolved with 
smoothed theoretical PSFs 9 to match the observed res- 
olutions: 6", 18" and 40" for 24, 70 and 160 /mi, re- 
spectively. All the fittings shown below are restricted 
to the data points between 100 to 700 AU where high 
signal-to-noise data are available. We first attempted to 
fit the radial profiles using a single grain size, but could 
not obtain a reasonable fit simultaneously for all three 
wavelengths. In fact, the [24]— [70] color indicates that 
the disk is composed of a hot component that emits ef- 
ficiently at 24 /j,m; while the comparable disk size at 70 
and 160 /im indicates that the disk is also composed of 
a cold component that emits efficiently at longer wave- 
lengths. 

We are able to find a reasonable fit for all three ra- 
dial profiles by assuming the disk has an inverse radial 
surface density distribution (p = 1), and is composed 
of grains with sizes: a± = 2 /im and ai = 18 /im, with 
the former size accounting for 98.5% of the disk particles 
by number. The fits are shown in Figure 1121 Assum- 
ing the disk extends to 1000 AU, the total mass of the 
material that emits at 24, 70 and 160 /im is 2.9xl0~ 3 

9 The observed PSFs are found to match well with smoothed 
theoretical STinyTim PSFs. If the object is well resolved compared 
to the instrument beam, convolving with a gaussian with a FWHM 
equal to the instrumental resolution is sufficient. However, if the 
object is not very well resolved, convolving with a gaussian will 
not match the observed data since MIPS PSFs are not gaussian- 
like especially at 160 fim. 
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Fig. 12. — The observed radial surface brightness profiles for 
Vega debris disk are plotted as open circles with colors representing 
different bands (24 fim: blue, 70 /an: green, and 160 fim: red). 
Model profiles from the large grain (a = 18 fim) component and 
the small grain (a = 2 /jm) component are plotted as dashed lines 
and dashed-dotted lines, respectively, with colors indicating the 
different bands. The sum of the two components is plotted as 
solid lines with colors representing the different bands. The 850 
fim emission from the two-component grain model is also plotted 
as a purple solid line for comparison to data from JCMT/SCUBA 
(purple open circles). 



M®, with 92% of the mass from the large grain compo- 
nent. However, the contributions of the small and large 
components to the total grain surface area are compara- 
ble. Note that the poor quality of the fit at the central 
part of the 160 /im profile is not sensitive to variation 
of the model parameters. The deficit at the central 160 
/jm model profile might suggest the existence of a unre- 
solved cold component close to the star. To further con- 
strain these parameters and estimate their uncertainties, 
we searched large regions of parameter space, comput- 
ing a x 2 statistic at each point. The search covered the 
following four-dimensional hypervolume: for the surface 
density power index, p, we searched from 0.5 to 3.0 with 
an interval of 0.1; for the number fraction of the small 
grains (fi), we searched from 10% to 99.9% with loga- 
rithmic spacing; for the grain size in the component 1 
(small grains, di) and component 2 (large grains, 02), 
we searched from 0.3 to 52 /im with logarithmic spac- 
ing. At 90% confidence, we conclude that p — 1.0 ± 0.2, 
fl = 98±£j %, ai = 2.0 ± 0.7 pm, and a 2 = 18tg 2 fim. 
The total mass estimated from the % 2 fitting is not well 
constrained (2.9±2° xlO -3 M e ) since the mass estimate 
is very sensitive to the specific size of the large compo- 
nent. 

We also fit the radial profiles using a grain size dis- 
tribution, n(a) oc a 9 , with minimum and maximum size 
cut-offs, a m i n and a max . The parameter space searched 
for the best fit encompassed values of q ranging from -1.0 
to -4.0, a m i n from 0.5 to 3.0 fim, and a max from 10 to 
60 fim. We fixed the surface number density index, p, 
at 1 because we found that the model profiles drop too 
quickly for p >1 in the previous two-component \ 2 fit- 
ting. At the 90% confidence level, we found a best fit with 
q = -3.0±0.5, a m in = 1-OtoJ M m , and a ma x = 46 ± 11 
fim. The best fit is shown in Figure El With these 
parameters, the total dust mass detected by MIPS is 
2.8±l.lxl0~ 3 M e . 

Our two-size and power-law distribution models both 
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Fig. 13. — Similar to Figure IT2l but using the model profiles 
from a grain size distribution: 11(a) oc a -3 with a minimum size 
cut-off a m i n = 1.0 fim and maximum size cut-of f am ax = 46 fim. 
The symbols and lines are the same as in Figure [H] 



fit the surface brightness profiles well at all three MIPS 
bands and can be used to predict the surface brightness 
profiles at 850 fim. The modeled 850 fim profile is shown 
as a purple line in Figures El an d El after convolution 
with a Gaussian beam-size of 16". For comparison, the 
reprocessed 1998 SCUBA archival data are also plot- 
ted as purple open circles after subtracting a ~5 mJy 
point source (representing the stellar photosphere) and 
smoothing with an additional Gaussian beam-size of 7". 
The total modeled flux emitted at 850 fim based on the 
two models is less than 3.5 mJy, much lo wer than what 
was f ound in the SCUBA data, 91±8 mJy ^Holland et alJ 
2005). 

To test the sensitivity of the model fits to the assump- 
tion of astronomical silicates, we also fit the radial pro- 
files using amorph ous carbon grains (p = 1.85 g/cm 3 , 
iZubko et aTl lT996). Similar to the silicate grains, a tem- 
perature distribution is computed for each of the grain 
sizes. Without changing the grain size distribution, a sat- 
isfactory fit is automatically obtained at all three MIPS 
wavelengths (Figure EJ- This is because the resultant 
temperature distributions in this environment are very 
similar between silicates and amorphous carbons (Figure 
111(1 . The total mass required to fit the MIPS data is 
1.5xl0 -3 M®, ~1.9 times lower than the silicate grains 
due to the reduced grain density in amorphous carbon. 
However, we rule out the possibility that the disk con- 
tains solely amorphous carbon as in that case the emis- 
sion at 850 fim would extend well beyond the observed 
SCUBA emission. 

A mixture of silicate and carbon grains is probably 
a more realistic dust grain model since both silicate and 
carbonaceous materials have been found in the Solar Sys- 
tem. Without spectral features to identify grain composi- 
tion, we simply assume that the particles in the disk con- 
sist of 70% silicates and 30% amorphous carbon grains. 
The best fit of this admixture model is displayed in Fig- 
ure !15l with a grain size distribution n(a) cx a - 30 ± - 6 an d 
a minimum size cut-off 3.2±0.8 pm and a maximum size 
cut-off 29±14 pm. The derived total mass in the disk is 
2.6±1.5xlO~ 3 M® with ~18% from the amorphous car- 
bon grains. The total emission at 850 pm is ~19 mJy 
assuming an aperture of 90". 

The independence of total mass on grain size distribu- 
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Fig. 14. — Same as FigureEDbut using amorphous carbon grains 
instead of silicates. At the MIPS wavelengths, the modeled profiles 
are very similar to the ones using silicate grains due to the simi- 
lar temperature distributions and emission efficiency. However, for 
a particle with a = 2 /im, the emission efficiency (Q a bs) a t 850 
(im for amorphous carbon grains is twice larger than the one for 
silicates. As a result, the modeled emission at 850 fim using amor- 
phous grains is much brighter and more extended compared to the 
emission using silicates. 



tion can be understood as a consequence of the drop in 
emission efficiency in proportion to the grain radius for 
grains much smaller than the wavelength; i.e., Q a bs °t a> 
when a <C ^ (Rayleigh limit). As a result of this ef- 
fect, the radiation goes as surface area times radius, or 
as a 3 . The number of grains required for a given surface 
brightness is therefore proportional to a~ 3 . The mass 
of a grain goes as volume, or a 3 . Therefore, the prod- 
uct of the required number of grains and the mass per 
grain is roughly independent of grain size. For grains 
with a <C A, the total mass of grains required to produce 
a given radiometric signature is roughly independent of 
a. Hence, the system mass derived from the MIPS data 
is insensitive to the specifics of the grain model. Since 
the temperature distributions for silicate and amorphous 
carbon grains are similar, hereafter we assume a pure sil- 
icate composition for our disk model, for simplicity. 

Our modeled flux at 850 is ~4-45 times lower 
than what was found in the SCUBA data. This implies 
that a third component is needed to account for the 850 
fim emission. We believe MIPS is seeing small grains 
blown out by radiation pressure while SCUBA is seeing 
larger grains that are gravitationally bound to the star 
and could be the grinding source of the small grains. The 
MIPS and SCUBA observations are then distinctly dif- 
ferent grain populations. To further test this scenario, 
we assume a ring structure in the disk, containing sili- 
cate grains with a = 215 /im. We do not try to model 
the ring's physical structure; instead we only try to find 
a structure where its resultant emission profile is consis- 
tent with the 850 /im profile and does not violate the 
24, 70 and 160 fim profiles. We achieved a fit with a 
ring structure that has a constant density from 86 to 100 
AU, and drops as r~ 2 from 100 AU to a cut-off radius 
of 200 AU. The radial profiles for this three-component 
model are shown in Figure ^j] In this model, the 24 
/tin profile is dominated by the small (2 /im) grain pop- 
ulation in the disk. The 70 /im profile, while dominated 
by emission from the large (18 /im) grains in the disk, 
has a significant contribution (<200 AU from the star) 



10.0000 




200 400 600 800 

Radius (AU) 

Fig. 15. — Same as Figures 1131 and \l4l but using a mixture of 
amorphous carbon (30%) and silicate grains (70%). The emission 
from amorphous carbon grains is plotted as dashed lines while the 
emission from silicate grains is plotted as dashed-dotted lines. The 
best fit grain size distribution is n(a) oc a -3 with a minimum size 
cut-off a m i n ~3.2 /im and maximum size cut-off a m ax ~29 /im. 



from the very large (215 /im) grains in the ring. Within 
this radius, the 215 /im grains contribute ^30% of the 
total emission. At 160 /im, the large grains in the disk 
and the very large grains in the ring contribute roughly 
equal amounts of emission. 

Having single-size grains in the ring is not mandatory. 
For silicate grains, we found that the modeled 24, 70 
and 160 /im radial profiles remain unchanged for grains 
with size larger than 180 /im since Q a bs ~1 for very 
large grains and their corresponding temperature struc- 
ture given the same heating source is similar (the canon- 
ical temperature distribution). Therefore, the total mass 
in the ring only changes by a factor of a 1 since the emis- 
sion goes as the grain surface area (a 2 ). Due to the uncer- 
tainty i n the absorption coefficient in the su bmillimeter 
regime l|Hildebrandlll983tlPoTlack et al]ll994D . we do not 
try to constrain the particle size distribution in the ring. 
The mass in the ring is at least a few tenths of a lunar 
mass (on order of 10~ 3 M®), con sistent with the previ- 
ous submillimeter measur ements ijZuckerman fe Becklinl 
lllMlHoiiand et alJll99i| . 

5.2. Stochastic Grain Heating 

The relatively constant [24]— [70] color temperature 
with radius suggests an alternative class of model in 
which tiny grains are heated stochastically by absorb- 
ing a single photon from the star. Stochastic heating 
occurs when the assumption that the temperature of a 
grain is determined by time-averages of the absorption 
and emission rates breaks down. This can occur when 
the the time interval between the absorption of photons 
with energies comparable to or larger than the grain heat 
capacity is much longer than the grain cooling timescale, 
a condition that can be realized in very small grains (< 
a few hundred A). To test the hypothesis of stochastic 
heating in the Vega system, we have computed t he com- 
plete time evolution of grain temperature (eg., iKriigel I 
2003) for small grains exposed to the radiation environ- 
ment of Vega. We have chosen a distance of 600 AU, 
characteristic of the debris disk extent, for the calcula- 
tion. 

For grains less than about ~50 A, the cooling timescale 



10 



Su et al. 



10.0000 




total emiss* 


Dn : 


1.0000 


" @- 




1 


0.1000 


^° 




-. 


0-0100 


:I I J 






0.0010 




\ y d (smoll + lorge) =2.0e-03 M® ( 
1\ M d (very large) =2.2e-03 M$ 


' — --^Cp ~ : 


0.0001 









200 400 600 800 

Radius (AU) 

Fig. 16. — Three component model fit with small (2 fim) and 
large (18 fim) grains in the disk, and very large grain (215 fi m) 
in the ring. The symbols and lines are the same as in Figure [Tzl 
with colors indicating wavelengths: 24 fim: blue, 70 fim: green, 
160 fim: red, and 850 fim: purple. The model emission from the 
215 fim grain in the ring is plotted as dotted-open-squares with 
color representing different wavelengths. The solid lines are the 
total model emission from the three-component model. 



was computed to be of order a few hundred seconds. 
For grains between 30-50 A, the absorption time scale 
varies from a few tens of seconds down to a fraction 
of a second in Vega's radiation field. With the heat- 
ing time scale much less than the cooling timescale, the 
temperature distributions are very strongly peaked near 
the equilibrium temperature and the emission from such 
grains can be very well described as an equilibrium pro- 
cess at a single temperature. Of course, any grain larger 
than 50 A will also satisfy the assumptions of equilibrium 
heating. At 20 A, the heating timescale has lengthened 
so that larger temperature fluctuations are observed and 
the temperature distribution of the grain has a significant 
hot tail. However, the emission from the grain is still 
very similar to that expected from an equilibrium pro- 
cess. In order to observe significant radiometric effects 
of stochastic heating in the Vega environment, grains of 
sizes <15 A are required. In this small regime, the heat- 
ing timescales are of order > 10 2 seconds and the grain 
will be stochastically heated. 

A large number ^10 A grains (which absorb ineffi- 
ciently in the optical and near UV) would be required 
to absorb sufficient energy from Vega to account for the 
infrared excess. As an illustration, we carry out a rough 
order-of-magnitude estimate. We assume that the grain 
absorption efficiency for a given wavelength of light goes 
as the geometric cross section times the ratio of grain ra- 
dius to wavelength. Then, typical stochastically heated 
10 A grains will have absorption efficiencies well below 
1% of the geometric cross sections (this estimate is con- 
firmed with real grain optical properties). If we put the 
grains at a typical distance of 600 AU, and assume a den- 
sity of 3 g/cm 3 , then we find that the mass of very small 
grains required to absorb the fractional luminosity of the 
star that is re-emitted in the infrared is within an order of 
magnitude of the mass required in grains in our preferred 
fit using large grains. It is not clear how to create the 
necessary numbers of such tiny grains without produc- 
ing substantial numbers of larger ones that would have a 
strong effec t on the radiometric properties of the system. 
In addition. lArtvmowiczl l)1988j) has calculated the effects 



of radiation pressure on such small grains around an A 
star and finds, because their absorption cross sections 
for stellar photons are very small, that they may not be 
ejected from the system by radiation pressure. Thus, 
the radial extent of the Vega system might be difficult to 
explain with tiny grains. The upper li mit to the polar- 
ization of the debris system obtained bv lMauron fe Dole! 
( 1998) also argues strongly against a system dominated 
by very small grains. It therefore is unlikely that a disk 
dominated by a population of small stochastically heated 
grains can explain all of the available Vega observations; 
nor is it necessary to invoke such a population as our 
preferred model fit using large grains is fully consistent 
with the observations. 

5.3. Summary of the Modeling 

A range of models using amorphous silicate and/or car- 
bon grains of sizes from ~1 to ~50 /im can fit the infrared 
radiometric behavior of the disk out to ^800 AU. The 
exact minimum and maximum grain size limits depend 
on the adopted grain composition. However, all these 
models require a r" 1 surface number density and a total 
mass of 3.0±1.5x 1O _3 M0 in grains. Models that also fit 
the submillimeter data additionally require an inner ring 
(near ~100 AU) with largr grains (>180 /im) and a total 
mass >10~ 3 M e . 

6. THE ORIGIN OF THE DEBRIS 

One might propose that the dust we see at 100-800 
AU is produced in situ from a highly extended Kuipcr 
Belt around Vega. Given that Vega is 2.5 times as mas- 
sive as the Sun, it is plausible that it could have a more 
extended population of such objects than does the solar 
system. Following simple mass scaling and assuming a 
similar surface density of large bodies in the disks, Vega's 
Kuiper belt would extend only about 1.6 times farther 
than the Sun's (45-55 AU). If Vega's Kuiper belt ex- 
tended to 500 AU, its mass would be over 100 times that 
of the Sun's Kuiper Belt (^IOOA/q). This mass, and the 
implied extent, seem implausible. Further, the distribu- 
tion of dust that would result from in situ production 
in such an extended Kuiper Belt almost certainly would 
not replicate the smooth surface brightness distribution 
we observe. We therefore disregard the hypothesis that 
the debris disk around Vega originates from an extended 
Kuiper Belt. 

The radial symmetry of the Vega disk strongly suggests 
the star system itself as the source of the dust. If the dust 
were of interstellar origin, it seems likely that it would 
be inhomogeneously dis t ribute d around the star. Also, 
lArtvmowicz fc Clampinl i|1997f ) show that the avoidance 
distance around Vega (the distance at which small grains 
would be repelled by radiation pressure), is ^3600 AU; 
therefore, small interstellar grains would not be present 
where we see the dust, inside about 1000 AU. Further, 
the radial density distribution of the dust is consistent 
with that expected from dust being ejected by radiation 
pressure. The debris we detect with MIPS is unlikely of 
interstellar origin. 

The large extent and the radial and azimuthal distri- 
butions of Vega's disk are all consistent with a model in 
which the dust originates at a distance ~100 AU from 
Vega, and is blown out to larger radii by radiation pres- 
sure. A further consequence of this model, as will be 
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shown below, is that the dust we are seeing in the Spitzer 
data originated in an event that took place in the rela- 
tively recent past. These conclusions do not apply to 
the larger particles in the ring structure around Vega, 
which are seen in the submillimeter but have much less 
contribution to the flux in the MIPS bands. 

The grain-sizes we deduce from our model fits to the 
data, ranging from ~1 to ~50 /im, are consistent with 
their being pushed out of the system by radiation pres- 
sure, if we assume the larger grains have significant 
porosity or are highly non-spherical. The ratio of the 
radiation force to the gravitational force on a grain is 
parameterized by (3 = i^Q^^p , where L„ is the stellar 
luminosity, <Q pr > is the radiation pressure efficiency of 
the grain averag ed over the stellar spectrum (< Q pr >~l 
when a »1 urn. lArtvmowicz fe Clampmlll997|) . G is the 
gravitati onal constant, M» is the stel lar mass, c is speed 
of light {Burns. Lamv. fc Soterlfl979l) . For (3 > 0.5, ra- 
diation pressure will drive the particle out of the sys- 
tem; for j3 < 0.5, the particle will spiral into the star 
due to Poynting-Robertson (P-R) drag. We adopt p — 
3.5 g/cm 3 for solid "astronomical silicate" grains. Us- 
ing L» = 60Lq and A/* = 2.5Af Q for Vega, we find 
02nm — 1-95 and ^lis^m — 0.22. The "small" grains 
in our model will be ejected, while the "large" grains 
will spiral in under the influence of P-R drag. In order 
to explain the presence of the larger (~ 20 /im) silicate 
grains at large distances from Vega, we require their den- 
sity to be reduced significantly. They may be aggregates 
of smaller grains, have significant porosity, or be highly 
non-spherical. It is also possible that the smaller (~ 
2^m) grains are under-dense, but that is not required to 
explain their presence at such large distances. 

Aggregate grains are seen in our own solar system (e.g., 
cometary dust and zodiacal dust. lCriin et aLll2001fl . and 
are a natural consequence o f growth through coagulation 
in th e early Solar Nebula (Kim ura. Okamoto. fc Mukail 
2002). Radiation pressure depends not only on the size 
and material of the dust grains, but also on their shapes 
and structures. Porosity certainly can decrease the bulk 
density p, therefore increasing the [3 value. In order for 
a e ff — 18 /im grains to have (3 > 0.5, their porosity has 
to be greater than 0.56 (or p e ff ~ 1.54 g/cm 3 , which 
is not unreasonable). If our large grains are very non- 
spherical (plates, needl es), they could also be ejected by 
radiation pressure (see lll'ln fc VoshchinnikovlflQM for 
example) . 

Particles escaping under the influence of radia- 
tion pressure are lost from the system on timescales 
of order the orbital period in the source region 
ijKrivov. Mann. fc Krivoval l2b00") . Considerations of en- 
ergy conservation lead to the conclusion that such par- 
ticles reach a terminal radial velocity given by v r ~ 

\j r GA t *i ^ ~ ^' wnere r imtiai is the distance where the 
grain is released or produced {Moro-Martin fc Malhotral 
|2005|) . This implies velocities for the escaping particles 
(assuming (3 =1) in Vega's disk in excess of 5 km/s, if the 
source region is around 100 AU, consistent with the loca- 
tion of the submillimeter dust ring. The residence time 
for the particles is roughly t residence ~ Rdisk/v r < 1000 
years for a disk radius of 1000 AU. Clearly, if the dust 
we are seeing in the Spitzer images is being blown out 



by radiation pressure, the individual dust grains are of 
recent origin. 

For our hypothetical very large grains (a >180 /im) in 
the inner ring at 86-200 AU that are seen in the submil- 
limeter data, we find that (3 ~0.02, implying very long 
residence times in that region. They therefore could be 
the source for the dust we see at much larger distances, or 
could be associated with a population of larger asteroidal 
bodies analogous to our own Kuiper Belt, which could be 
the ultimate source of the dust at large distances. Cas- 
cades of collisions starting with encounters between these 
larger bodies in this ring could produce the small grains 
that are blown out by radiation pressure force. 

The short residence time for small grains (a <50 ptva), 
combined with the dust mass we derived in section 5 (see 
Figures H 31 16|) . places a lower limit on the dust produc- 
tion rate in the Vega system. Taking the dust mass in 
small grains to be 3.0 x 10 -3 M® and a residence time of 
100 years, the require d dust production rate is ~6x 10 14 
g/s. iLisse et "all l)1999f) estimate an overall dust produc- 
tion rate of 2xl0 8 g/s for the Comet Hale-Bopp. If the 
dust is generated by cometary activity, the Vega debris 
disk requires 3 million Hale-Bopp- like comets to account 
for the dust production rate, which is impossibly high. 
If we further assume that dust loss and production have 
been in approximate balance over the life of the Vega sys- 
tem (350 My), we find that the mass of the dust lost over 
that time is ~ 1100 M® = 3.3 Mj. This is a lower limit 
as we have assumed the longest plausible residence time, 
and limited the calculation to only those small grains we 
are sensitive to; presumably there are other grains lost to 
P-R drag in the inner part of the disk. Further, to pro- 
duce 3.3 Mj worth of dust, a much larger mass of parent 
bodies would have to be present in the disk. Assuming a 
power law size dist ribution propor tional to a -3 ' 5 , where 
a is the radius (see Dohnanvil ll969ri . the mass of the par- 
ent bodies would be two to three orders of magnitude 
greater, implying a total mass of 330-3300 Mj (or 0.3-3 
Mq) for Vega's debris disk. 

It seems very implausible to have an initial inventory of 
330-3300 Mj of dust when Vega was born. Only a frac- 
tion of the initial disk mass is presumably incorporated 
into asteroids; most of it could have been sequestered 
in giant planets, accreted onto the star, or ejected from 
the system. Besides, debris disks do not form as solid 
material only; they start out as massive protoplanetary 
disks, with 100 times as much gas mass as dust. The 
disk masses of Herbig Ae/Be stars (what Vega would ap- 
pear as, when first formed) are on the order of a few 
tenths of a M^, measured in the millimeter contin uum 
{Hillenbrand et alJll99l IMannings fc Sargentill997j) . It 
is very unlikely that this high dust production rate is in 
a steady state throughout Vega's life. 

The above arguments suggest that the disk we imaged 
with Spitzer is ephemeral. Such a disk could be pro- 
duced as the result of the disruption of a large aster- 
oidal or cometary body and subsequent collisional cas- 
cade as pieces of the disrupted body collide with other 
bodies in the source region. The dust mass we observe, 
3.0xlO~ 3 Af e , would form an object ~1000 km in radius 
(the size of Pluto) if collected together. Such a large 
object is unlikely to unde rgo a disruptive collision {e.g. 
Kenvon & Bromlev 2004). We hypothesize instead that 
a more modest sized object was disrupted, and the ensu- 
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ing collisional cascade, interacting with the background 
population of asteroidal objects, produced the observed 
dust. Small bodies were probably the ultimate source of 
most of the d ust mass, although regolith ej ection from 
larger bodies (Kortcnk amp fc Dermottlll998fl could also 
have contributed significantly. Thus we favor a scenario 
in which the dust we observe is young because it is be- 
ing blown out of the system by radiation pressure, but 
was produced over some finite interval of time, consistent 
with a collisional cascade subsequent to a large disruptive 
co llision. This scenar io is also consistent with the results 
of iRieke et all l)2005t) . who find that episodic collisions, 
and not steady-state production, best explain dust seen 
around a large sample of A stars. 

7. CONCLUSION 

We have presented images of the Vega debris disk at 
24, 70 and 160 /im obtained by MIPS. The disk is well 
resolved with 1-er detection radii of 43", 70", and 105" at 
24, 70 and 160 /im, respectively. The disk appears cir- 
cular and smooth at all three wavelengths. No clumpy 
structure was found. The surface brightness of the disk 
follows a simple r~ 3 (inner disk) or r~ 4 (outer disk) 
power law at all three wavelengths, implying the disk 
density structure is simple and smooth. A region lack- 
ing material that emits efficiently at MIPS wavelengths, 
was inferred at a radius of 11"±2" from a radial profile 
analysis. 

We investigated the nature of the Vega disk with model 
fitting. We found that an axially symmetric, face-on 
disk with r^ 1 surface number density can explain the 
observed radial profiles at all three wavelengths simulta- 
neously. The observed radial profiles can be well repro- 
duced with a two-component grain model consisting of a 
mixture of small (a = 1.3-2.7 /im) and large (a = 12-30 
jum) silicate grains with the former size accounting for 
^98% of the particles in number. The observed radial 
profile can also be reproduced with a power-law grain 
size distribution, n(a) ~ a ~ 3 - 0±0 - 5 with a minimum size 
cut-off 0.7-2.3 /im and maximum size cut-off 35-57 /im. 
The minimum and maximum grain sizes in the disk are 
constrained by the requirement that the model simulta- 
neously reproduce the 24, 70, and 160 fim radial profiles. 
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Silicate and carbonaceous grains have very similar ra- 
dial temperature profiles in the radiation environment 
around Vega. Therefore, while the size distribution is 
well constrained by the surface brightness profiles, the 
grain composition is not. While the composition of the 
dust grains in the disk is unknown, it is fully consistent 
with a mixture of silicate and carbonaceous dust. Re- 
gardless of the exact composition, the total mass of the 
dust seen by MIPS is ~3.0xl0 -3 M ffi assuming an outer 
radius of 1000 AU. 

The total flux at 850 /im emitted by the grains required 
to fit the MIPS data is much less than what was measured 
in the SCUBA data. Another dust component is needed 
to account for the emission at 850 /im. Emission from a 
ring composed of grains larger than 180 /urn can account 
for the observed 850 /im profile, while their contribution 
in the MIPS wavelengths does not change the disk model 
profiles. 

The ring-like structure in our disk model represents a 
reservoir for gravitationally bound planetesimals, where 
collisions can occur and generate small dust grains as 
debris. The dust production rate implied by our obser- 
vations, ~10 15 g/s, would require this asteroidal reservoir 
to be improbably massive, were the disk we observe with 
Spitzer in a steady state. We conclude that the disk is 
ephemeral, the consequence of a large and relatively re- 
cent collisional event, and subsequent collisional cascade. 
Radiation pressure sets the small dust grains produced in 
these collisions on unbound hyperbolic trajectories and 
they stream out through the disk and leave the system 
with a time scale on the order of ~1000 yr. The large 
extent of the disk as seen by MIPS is consistent with the 
dust being ejected by radiation pressure. 
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TABLE 1 
Observing Log 



Array 


Field Size 


Pixel Size 


Exposure 


# of Pos. 


Eff. Exp. 


1-a per pixel 






arcsec 


exp x eye. 




sec 


mjy/arcsec 2 


A24 


small 


2.55 


3 s x 1 


4 


120~126 a 


0.01 


A70F 


small 


5.20 


3 s x 1 


12 


132~151 


0.13 


A70 


large 


9.98 


3 s x 8 


1 


195~255 


0.04 


A160 


large 


16x18 


3 s x 4 


4 


38~66 


0.02 



a The typical effect exposure time is ~ 63 sec for the central 6" region due to the 
saturation. 
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